Journal of Nonlinear Mathematical Physics 2001, V.8, N 1, p^ioi 



Article 



Integrability and Explicit Solutions 

in Some Bianchi Cosmological Dynamical Systems 

J CHAVARRIGA 1 and I A GARCIA 2 

Departament de Matemdtica, Universitat de Lleida 
Avda. Jaume II, 69, 25001, Lleida, Spain 

1 E-mail: chava@eup.udl.es 

2 E-mail: garcia@eup.udl.es 

T— I ! 

Received April 6, 2000; Revised August 2, 2000; Accepted September 18, 2000 

(N : 

Abstract 

« , 

^ ^ ■ The Einstein field equations for several cosmological models reduce to polynomial 

systems of ordinary differential equations. In this paper we shall concentrate our 
attention to the spatially homogeneous diagonal G2 cosmologies. By using Darboux's 
I I theory in order to study ordinary differential equations in the complex projective 

• plane CP 2 we solve the Bianchi V models totally. Moreover, we carry out a study of 

t**"* ' Bianchi VI models and first integrals are given in particular cases. 

<n : 
t— 1 . 

^ ! 1 Introduction 

The Einstein field equations for a number of classes of cosmological models have been 
*q ' written as polynomial systems of ordinary differential equations. In recent investigations 

some authors look for new exact solutions of these differential systems due to their cosmo- 
?h . logical interest. Most of these solutions araise either due to the existence of a symmetry, 

or under a restriction on the parameters of the cosmological model. Homogeneous cos- 
. ^ ' mological models have important applications in the theory of evolution of the universe. 

^ ■ In particular they can be applied to the theory of explosions of stars, formation of gala- 

xies, etc. The equations of general relativity are reduced to systems of a finite (but quite 
large) number of ordinary differential equations. In the last two decades these dynamical 
systems were and still are being analyzed by means of traditionals analytic and numerical 
methods. However the dynamical systems in question are so complex and the development 
of effective methods for their analysis becomes especially urgent. 

In this paper we shall concentrate our attention to the diagonal G2 cosmologies, i.e. 
those spacetimes which admit an Abelian G% comprising two hypersurface-orthogonal 
Killing vector fields (see for instance Q). In fact we only study spatially homogeneous 
diagonal G2 cosmologies. It is assumed that the Einstein field equations are satisfied and 
that spacetime admits a perfect fluid with linear equation of state 

2 

p=(<y-l)p, -< 7 <2, (1) 
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where 7 is an arbitrary parameter in the specified range of values which contains the 
physically significant values 7 = 1 (corresponding to dust) and 7 = 4/3 (for which the 
fluid is interpreted as radiation). In all the paper we use geometrized units with c = 1 
and 8ttG = 1 and the sign convention of MacCallum B . 

In Section 2 we describe the mathematical techniques used for the analysis of polynomial 
differential equations that we deal in this work. A summary is given about the Darboux's 
theory in order to study ordinary differential equations in the complex projective plane CP 2 
and, in short, its utility for the analysis of the socalled degenerate infinity systems. This 
name is due to the fact that in the Poincare compactification of these kind of differential 
systems all the equator of the sphere S 2 , i.e. the infinity of the affine plane, is filled up 
of critical points, see Sotomayor In Section 3 we show the equations that govern the 
spatially homogeneous diagonal G2 cosmologies and the Bianchi models. In Section 4 we 
solve the Bianchi V models totally and, finally in the fifth and sixth sections we carry 
out a study of the Bianchi VI cosmological models. First integrals are given in particular 
cases and a canonical form of the equations allows us to obtain the formal solution. 



2 Polynomial foliations in CP 2 

In this section we consider the first order differential equation 

Q{x,y)dx-P{x,y)dy = 0, (2) 

where P and Q are polynomials which in the following we can suppose with coefficients in 
the complex field C. We say that m = maxjdeg P, deg Q} is the degree of polynomial dif- 
ferential equation (||). When we consider it necessary, we introduce a "time" parameter t, 
which enables us to express (Q) as a first order autonomous system in the affine plane 

X = ~d~t = P<yX, V ^ V = ~c(t = ®( X ' ^' ^ 

Following Darboux ||, we introduce the homogeneous coordinates (X,Y,Z) of the 
complex projective plane CP 2 such that 

X Y ^ 

x= z' y= z' (4) 

which allow us to express differential equation (|2|) in the form 

L(X,Y,Z)(Y dZ - ZdY) +M(X,Y,Z)(Z ,dX-XdZ) = 0, (5) 

where L(X,Y,Z) = Z m P(X/Z,Y/Z) and M(X,Y,Z) = Z m Q(X/Z,Y/Z) are homoge- 
neous polynomials of degree m. Note that Z = is an invariant straight line of this 
equation. 

More generally (see Q), we can consider a differential equation defined in CP 2 as an 
equation of the form 



L(Y dZ-Z dY) + M(Z dX-X dZ) + N(X dY -Y dX) = 0, 



(6) 
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where L(X, Y, Z), M(X,Y, Z) and N{X,Y, Z) are homogeneous polynomials such that 
deg L = deg M = deg N = m. Introducing the auxiliary independent variable T, we can 
also interpret equation (g) like a dynamical system 

X' = ^ = L(X,Y,Z), Y ' = ^ = M(X,Y,Z), Z> = ^ = N(X,Y, Z), (7) 

over C 3 \{0}. Such homo geneous system induces on CP 2 the polynomial foliation given 
by (|6|). In fact, differential equation (|6|) has not a unique expression. 

Remark 1. If we take L = L- AX, M = M — AY and N = N — AZ instead of L, M 
and N respectively into equation (||), where A(X,Y,Z) is a homogeneous polynomial of 
degree m — 1, then differential equation Q remains invariant. 

Remark 2. One of the most obvious questions to ask is which is the relationship between 
the parameters t and T of the affine differential system (|3|) and the projective differential 
system (0) with N = 0. A straightforward differentiation shows that 

^ = Z m - 1 [L(x, y, I) - xN(x, y, 1)] , ^ = Z m ~ l [M(x, y, I) - yN(x, y, 1)] . 
On the other hand, by the chain's rule we have 

^ =P(xv )± ^L = Q(xy) — 
Finally, identifying the last two previous expressions when N = we conclude that 

m 

A polynomial differential equation (|2|), where P(x,y) = Pk{ x ,y) an d Q(x,y) = 

k=0 

m 

Yl Qk( x , y) is the developement in homogeneous components of the polynomials P and Q 

k=0 

is said to be degenerate infinity if it verifies the condition xQ m (x,y) — yP m (x,y) = or 
equivalently P m (x,y) = xA(x,y) and Q m (x,y) = yA(x,y) for some homogeneous polyno- 
mial A{x,y) of degree m — 1. In order to clarify this terminology, we write equation (|2|) 
in the complex projective plane as equation @, where 

m— 1 



L(X, Y,Z) = Y J Pk(X, Y)Z m ~ k + A(X, Y)X, 



k=0 
m—1 

M(X, Y,Z) = Y, Qk(X, Y)Z m ~ k + ApT, Y)Y, 

k=0 

N(X, Y, Z) = 0. 

Now, from Remark 1, we can take L = L — AX, M = M — AY and N = N — AZ obtaining 
L(Y dZ -Z dY) + M(Z dX - X dZ) + N(X dY -Y dX) = 0, (9) 
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and being 



m— 1 

L(X, Y,Z)=ZJ2 P k(X, Y)Z m ~ k -\ 

k=0 
m— 1 

M(X,Y,Z) = ZJ2Qk(X,Y)Z m - k -\ 

k=0 

N(X,Y,Z) = -ZA. 

From the above equation we conclude that the line of infinity Z = is obviously an 
invariant straight line of (Jsj) and moreover this line is filled up in a dense sense of critical 
points since L(X, Y, 0) = M(X, Y, 0) = N{X, Y, 0) = 0. 

At this point we want to comment that it is more convenient, in order to obtain their 
critical points, write equation (^) into the form 

VdX + QdY + KdZ = 0, (10) 

being V = MZ - NY, Q = NX - LZ and 1Z = LY — MX. We say that (X , Y , Z ) G 
CP 2 is a critical point of equation (10) if and only if V(Xq, Yq, Zq) = Q(Xq,Yq, Zq) = 
TZ(X ,Y ,Z ) = 0. 

Finally, we recall that the complex projective plane CP 2 is compact because it is ob- 
tained by adding a line at infinity to the plane C 2 . Hence, when we extend the foliations 
from the affine plane to the projective plane and we study their integrability, this study 
is global in the above sense. 



3 Spatially homogeneous diagonal G2 cosmologies 

Homogeneous cosmological models are defined to be space-time manifolds thogether with 
a metric ds 2 satisfying Einstein's equations and invariant under some three-dimensional 
Lie group of transformations of the manifold. In cosmology the most important models 
are spatially homogeneous models for which the orbits of the Lie group are space-like 
(restrictions of the metric ds 2 to the orbits of the group are negative definite). 

Let S + and S_ be dimensionless shear components. In addition, the variables A and -/V x 
are dimensionless measures of the isotropic and anisotropic curvature of the hypersurfaces 
which lie orthogonal to the fluid flow. It is easy to show (see ||] and ||] for example) that 
the Einstein field equations for the spatially homogeneous diagonal G2 cosmologies reduce 
to the following system of polynomial ordinary differential equations 

E+ = (g-2)£ + -2iV 2 , E_ = (q - 2)E_ + 2AN X , 
N x =(q + 2S + )iV x , A = (q + 2£ + ) A, 

where the over dot denotes d/dr, i.e., derivative with respect to a dimensionlees measure 
of time r along the fluid flow lines (see PH| ). Moreover, the deceleration parameter q is 
given by 



^S^ + S 2 ) + -(3 7 -2)0, (12) 
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and the density parameter adopts the form 

tt = 1 - T? + - E 2 _ - A 2 - Nl . (13) 



For this cosmological model it is known that the Einstein field equations (11) admit 
the particular solution 

A£+ + iV x £_ = 0. (14) 

The constraint equation ( |i~4"l) allows the following possibilities (each one of them gives 
rise to a Bianchi model). 

• If A = N x = we have the model B(I). 

• Imposing S + = N x = the resultant model is B(V). 



• From the last two equations of (11) it is obvious that if N x ^ 0, then system ( |Tl|) 
admits the first integral A = kN x being k a constant. Therefore, constraint (|i~4]) 
implies £_ = —kT, + and the cosmological model is called B(VI). 

The Bianchi I models B(I) are expressed by the following system 

£+ = (<7-2)£+, £_ = (g-2)£_, (15) 



where q is given by ( |l2|) and Q = 1 — — S 2 .. These models are easily integrated to yield 
either £_ = kT, + , where k is an arbitrary constant or S + = 0. Therefore, in the phase 
portrait we only have singular points and straight lines. 

In |H the author carries out a qualitative study of all these Bianchi type models. On 
the other hand, Q examine the different Bianchi types and determine all of the linear and 
quadratic algebraic invariant curves which lie in the physical region of phase space. 



4 The solution of Bianchi V models 

The Bianchi V models B(V) are given by the following cubic polynomial systems 

£_ = (g-2)£_, A = qA, (16) 
where the deceleration q and the density Q, are 

q = 2Y?_ + -(37 - 2)0, O = 1 - A 2 - Y?__. (17) 

These models have been studied by many authors. In particular, it has been pointed 
out (see [||, 0] and |J) that the Einstein field equations reduce to elliptic integrals for 
certain values of 7. Moreover, in § the following time invariant for £_ 7^ is given 

A = Q S_e 2r , 

where Qq is an integration constant. 

Performing the change of variables £_ = y/x, A = ^fy we can reduce the degree of 
system ( |l6| ) obtaining the following quadratic system with degenerate infinity 



x = 3(7 - 2)x + xA(x, y), y = (37 - 2)y + yA(x, y), 



(18) 
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where A(x, y) = 3(2 — 7 )x + (2 — 3 7 )y. Applying the arguments of second section, we write 
system ( |l~8| ) in the complex projective plane as 

L(Y dZ - Z dY) + M(Z dX - X dZ) + N(X dY -Y dX) = 0, (19) 

where, making the rescaling of time dr = Z dT (see Remark 2), we have 

L = 3(7 - 2)X, M = (3 7 - 2)Y, N = -A(X, Y). 



From the first two equations L and M we obtain for equation (1£) the following first 
integral 

Hx(X,Y,Z) = x^- 2 Y^ 2 ~^\ 

which is a fourth degree homogeneous function. Moreover, as L, M and A~ are linear 
homogeneous polynomials which only depend on the variables X and Y, there exist con- 
stants ci, C2 and C3 not all zero such that c\L + c<iM + C3N = 0. Hence, we have for 
equation ( |l9|) the first integral H^X, Y, Z) = c\X + c^Y + c^Z. In short we have for 
instance the linear first integral H2(X,Y, Z) = X + Y — Z . Now, we consider for equa- 
tion (ED the zero degree homogeneous first integral defined like H\ (X, Y, Z) / (X, Y, Z) . 
Making the transformation (X, Y, Z) — > (x, y, 1) in the above first integral we turn to the 
; plane obi 

H(x,y) :-- 



afflne plane obtaining the first integral for system (18) 
fTi(a,y,l) _ x 37-2 2/ 3(2- 7 ) 
Hl(x,y,l) ~ {x + y-iy 
Finally, undoing the change of variables S_ = y/x, A = y/y we have for Bianchi V 



models (16) the following Darboux first integral 



^37-2^3(2-7) 

H(E-,A) = — 5 . (20) 

(S^ + A 2 -l) 2 

Moreover, from the equations of L and M we obtain the time parametrization 

X(t) = AT e 3(7 - 2)r , Y(t) = y e (37 - 2)r , 

where Xq and Y$ are real parameters. In addition, since A^ = —A(X,Y), we can integrate 
this equation and write the time parametrization of the variable Z as follows 

Zir) = - J A(X(t), Y(t)) dr = - J [3(2 - 7 )X(r) + (2 - 3 7 )F(r)] dr 

= - J [3(2 - 7 )Xo e 3 ^- 2 ) r + (2 - 3 7 )Y e( 3 ^ 2 ) r ] dr = X(r) + Y(r) + C, 

being C a real parameter also. Finally, we turn to the afflne plane M 2 by means of x(t) = 
X(t)/Z(t), y(r) = Y(t)/Z(t), and we have the general explicit solution of Bianchi V 
models (|l8|) given by S_(r) = ^/x(t), A(t) = \fy{r). In a more explicit way, taking the 
initial values E_(0) := S° = ^A: /(A:o + Y + C) and A(0) := A = y/Y / (X + Y + C), 
we have 




exp(3( 7 - 2)r) 

:p((3 7 -2)r)-l]+l' 

(21) 



SO f [exp(3( 7 - 2)r) - 1] + (^°) 2 [exp((3 7 - 2)r) - 1] + 1 



exp((3 7 - 2)r) 



(S°) 2 [exp(3( 7 - 2)r) - 1] + (,4 ) 2 [exp((37 - 2)r) - 1] + l" 
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5 Bianchi VI models 

The Bianchi VI models B(VI) are given by the following cubic polynomial system 

t + = {q-2)Y, + -2N 2 x , N x = (g + 2£ + )i\T x , (22) 
where the deceleration q and the density £1 are 

q = 2(l + k 2 )E 2 + + ^ -2)n, n= 1- (1 + k 2 ) (N$ + £+) . (23) 

These models admit the invariant algebraic curves N x = (Bianchi I solutions) and f2 = 
(Ellis-MacCallum vacuum solutions, see J/J p. 136). 
In H Hewitt shows that in the particular case 

, ( 8 



1 + ■ (24> 



system ( |22| ) possesses another invariant algebraic curve given by 
/ 3 7 + 2\ 2 Nl /3 7 + 2\ 2 



Furthermore, for this value of the parameter k, he shows that system ( p2j ) has the time 
invariant 



^(2-7)437+2^-8 = Qoexp j^_9 (2 _ 7)2 \ 



(26) 



But, with these three invariant algebraic curves, we can be able of construct for 
Bianchi VI models (|22] ) with restriction (|24| ) the integrating factor 

37+2 97-2 

i?(£+,iV x ) := iV~ 3 A 3(7-2)^6(2-7). (27) 

There exists another particular case of Bianchi VI model given by the constraint 
, ( 2 



1 + * = [*f=Z)> < 28 > 



for which system ( p2j ) admits the following new quadratic invariant algebraic curve 

r := (2 - 3 7 )(3 7 - 4) 2 + 4(4 - 3 7 )(3 7 - 2)£+ + 4(2 - 3 7 )S 2 h + 12(2 - 7 )iV x = 0.(29) 



This invariant curve allows us to obtain for Bianchi VI models ( p2| ) with restriction (|2 
an integrating factor given by 

2(4-37) 

R(E+,N X ) -.= n x ^- 2 r 3 ^- 43 ^- 2 fr 2 . (30) 

Thus, the first integral H associated to the integrating factor R is given by 
H(Z + ,N X ) = J [(q-2)Z + -2N 2 x }R(Z + ,N x )dN x +F(X + ), 
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where the function F is calculated as usual by the condition 
B M 

-R(X + ,N x )(q + 2X + )N x . 



3£_ 



Another point of view in order to study Bianchi VI cosmological models it is the 
transformation of system (22) to the projective differential equation (|J) by means of the 
change to homogeneous coordinates S + = X/Z, N x = Y/Z. Thus, after the application 
of Remark 1 we obtain that system (^) adopts the form 

L(Y dZ - Z dY) + M(Z dX - X dZ) + N(X dY -Y dX) = 0, (31) 

where 

L(X, Y, Z) = X -Z (3( 7 - 2)XZ - AY 2 ) , 

M(X, Y, Z) = l -ZY {AX + (37 - 2)Z) , 
N(X, Y, Z) = -ZA(X, Y, Z), 



and being A(X, Y, Z) := \ (l + k 2 ) [3(2 - 7) A 2 + (2 - 37)y 2 ] . Differential equation © 
admits the first integral TC(X,Y, Z) := Z 2 Q(X/Z,Y/Z) = h. So, taking into account this 
first integral we can reduce the dimension of equation (|3l]) and we obtain the following 
differential equation 

dX _ Ah + 4 (1 + k 2 ) X 2 + 3 (1 + k 2 ) ( 7 - 2)XZ - AZ 2 

~dZ ~ (1 + k 2 ) [(2 - 3 7 )/i -4(1 + k 2 ) X 2 + (3 7 - 2)Z 2 ] ' (32) 
which has the particular solution J-(X, Z) := —h — (l + k 2 ) X 2 + Z 2 = 0. The study of 



differential equation ( |32| ) allows us to obtain other integrable Bianchi VI models. Thus, 
for instance, if relation 



2 

„2 



1 + k -{^) • (33> 

is verified, then the straight line F\(X, Z) := AX + (37 — 2)Z = is an algebraic invariant 
curve of equation (|3^). Moreover, turning to the affine plane we obtain that Fi(E + , 1) = 
is an invariant straight line of the Bianchi VI models ( p2| ) with constraint (^). Now, with 
this new particular solution we can construct the integrating factor 

2(4-37) 97 _ 2 

2?(E+,JV X ) := N^- 2 (34) 

The first integral associated to the above integrating factor R is given directly by f(q + 
2E + )iv x i?(£_|_, iVx) A more detailed study of this first integral shows that there 

exists another quadratic invariant algebraic curve F% which is functionally independent 
of $7. In short, we have that 

^ 2 (£ + , iVx ) := (37 " 2) 4 (37 + 2) S + + (2 + 3 7 )^ + (3 7 - 2)iV 2 = 0, 

is a new invariant algebraic curve for Bianchi VI models with constraint (|33|) which allows 
us to obtain the following Darboux first integral 

H(Z + ,N X ) = Nl {l ~ 2) Vt 2+ ^'Fl {2 ^\ (35) 
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6 Canonical form of Bianchi VI models 

If we rescale the variable X by X = X\J\ + k 2 and time by dt = 3(7 — 2) (l + k 2 ) dT, we 
write differential equation (|32|) in the following canonical form 



X' = XZ + \ X T{X,Z), 

(36) 

Z' = X 2 + \ 2 F(X,Z), 

where the prime denotes derivative respect to t. Now J-(X, Z) := —h — X 2 + Z 2 is an 
algebraic particular solution of the canonical system (|36| ) and the parameters Ai and A2 
are given by 

A,= - 4 , M 



l(7-I)Vlt?' 3(7-2)' 

The integrable Bianchi VI models which arises from parameter restrictions (|24|), (|2£ 
and p^ ) correspond to integrable cases of system ( p6| ) with Ai + A2 = 1/2, Ai + A2 = — 1 
and Ai + A2 = respectively. 

For Ai + A2 = 1/2, the canonical system ( |36"| ) has two new invariant algebraic solutions, 
namely T\{X, Z) := \fh + X + Z and T%{X % Z) := \fh — X — Z. Now, performing the 
change of variable U = X + Z the canonical system (|36|) becomes 



X' = X{U - X) + Ai {-h + U 2 - 2UX) , 

1 (37) 
U> = \{U*-h). 



Therefore, by integration of the second equation of (37), we have that the general solution 
of the U variable is given by the parametrization 

U(t) = ^-eMVHt-to)]^ 
1 + exp[Vh{t - t )} 

where to is an integration constant. Moreover, the first equation of system (|37| ) is of Riccati 
type and in addition it has a particular solution X\(t) which we known from the invariant 
algebraic curve J r (X 1 ,U - X x ) = 0. In short, we have X x {t) = [-h + U 2 (t)] /2U(t) or 
more explicitly 

X\(t) = V/icosech[v / /i(t — to)], 

and hence we can obtain the parametrization X (t). Moreover, undoing the changes of 
variables and rescaling of the independent variable we obtain the parametrization X(T), 
Y(T) and Z(T) of the solutions in the complex projective plane in terms of the dilogarithm 
special function Li2(z) := f 2 ° t _1 log(l— t) dt. But the turn to the affine plane is so complex 
because the relationship among the physical time r and the projective time T is expressed 
as dr = Z(T) dT and this integration is difficult. 
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